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Abstract 

We introduce coarse-grained hydrodynamic equations of motion for diffusion- 
annihilation system with a power-law long-range interaction. By taking into 
account fluctuations of the conserved order parameter - charge density - we 
derive an analytically solvable approximation for the nonconserved order pa- 
rameter - total particle density. Asymptotic solutions are obtained for the case 
of random Gaussian initial conditions and for system dimensionality d > 2. 
Large-t, intermediate-t and small-t asymptotics were calculated and compared 

with existing scaling theories, exact results and simulation data. 
PACS: 05.60.+W, 82.40.-g, 61.30Jf 
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I. INTRODUCTION 



In recent years, the annihilation-diffusion problem has generated significant interest, 
both theoretical and experimental. The annihilation-diffusion problem usually corresponds 



most physical systems that one is interested in, in addition to the thermal diffusion and ki- 
netic annihilation, the particles A and B are charged and interact via a power-law long-range 
interaction (LRI), usually of Coulomb type, although here we will consider a more general 
LRI. Such physically important interaction can clearly strongly inffuence the annihilation 
dynamics by introducing an additional time scale in the annihilation process, and can lead 
to a new mechanism for slow dynamics. The two-species annihilation reaction with LRI can 
be studied in a variety of experiments. A thermal quench of a freely-suspended liquid crystal 
film from smectic-A to smectic-C phase is one experimental system where this annihila- 
tion process has been studied in great detail. In such experiments, immediately after quench, 
the singularities of the smectic director (2d vector) field appear as positive and negative vor- 
tices interacting (due to elastic forces) via a logarithmic potential. As time elapses, vortices 
of opposite sign slowly annihilate, exhibiting complex dynamics that is clearly and strongly 
influenced by temperature, initial particle distribution, and LRI that leads to an attraction 
between annihilating partners. Similar annihilation problems also appear in turbulent flow, 
superconductivity, spinodal decomposition, and many other condensed matter systems. The 
annihilation process is also relevant to coarsening of topological defects produced by sym- 
metry breaking field in particle physics models, after an early temperature quench due to 
the fast initial universe expansion, a process that is thought in part to determine the large 
scale structure of today's universe. 

It is well known that in classical chemical kinetics, density decay for both one-species 
and two-species annihilation is described by the kinetic rate-equation (see, e.g., 0] ): 



to that of the kinetics of particle density decay in the annihilation reaction A + A 







(one-species annihilation) and A + B 



(two-species annihilation). In the latter case, in 




(1) 
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with the large-t asymptotics given by: 

Pit) ^ (JCt)-' , (2) 

where /C is a reaction constant. Equation (|1]) completely neglects all spatial fluctuations 
and correlations of the particle density. Although the dynamics described by equations (|Tp 
- (0) is correct above an upper critical space dimension due, for d < due, both diffusion and 
fluctuations play an important role, substantially slowing down the density decay. It was 
indeed shown P], that for the one-species annihilation without LRI, the due = 2. That is, 
for systems of dimensionality less than 2, the density decay is given by: 

Pit) ^ poiDp^^/'k)-'/', (3) 

where D is diffusion constant and po is the initial particle density. Equation (j^) can be 
understood either by invoking a single length scaling arguments or by using a more 
elaborate Smoluchowski approach For the two-species case, it was shown (Refs. 0, 

that the upper critical dimension due = 4, and that for d < 4, the large-t asymptotics 



is: 



Pit) ^ ipo)'/'iDt)-'/\ (4) 

The decay law (P was conflrmed in several numerical simulations of one-, two-, and three- 
dimensional systems It is important to emphasize that in order to observe such 
power-law decay, it is necessary initially to have an equal number of positive and negative 
charges, distributed at random. If initial numbers of positive and negative charges are 
different, one should observe an exponential density decay to the non- vacuum equilibrium 
(see, e.g., Sokolov et al. Jl^)- If the system is well-stirred, i.e., long- wavelength fluctuations 
of charge density are suppressed, the decay law would also be different. We will not consider 
such cases in this paper. 

To better study the role of correlations in reactions without LRI, description in terms of 
secondary quantization operators of creation and annihilation was proposed. ||l5|Jl6[| Further 



developing this approach, Pehti proposed an "exact" (vahd to all orders in perturba- 



tion expansion) renormalization-group theory for the one-species annihilation, and Lee and 
Cardy ||18| suggested a similar approach to the two-species annihilation, both based on the 
rigorous master equation converted into a field-theoretic formulation. In the latter case, it 
was shown that the decay law (^) is an asymptotic one for large t and small po- It was also 
rigorously proven that the assumptions used in the derivation of equation (|) are correct 
only for d > 2, while for d < 2, a more elaborate renormalization group procedure which 
includes an effective noise with nontrivial correlations should be carried out. 

The addition of the LRI further complicates the picture. Equations of motion be- 
come strongly coupled via the nonlocal interaction, and only approximate solutions of 
these equations can be found. Up to date, only Coulombic |]19| systems in two dimensions 



{d = 2,n = 1, the case corresponding to particles interacting via a logarithmic potential 
on a two-dimensional substrate or film) were studied numerically |[r3|j20| -p2| , and results 



appear to be inconclusive. In these simulations, particle density exhibited a power-law de- 
cay p(t) ~ t-" with the exponent u varying from 0.79 ± 0.04 to 0.85 ± 0.05 The 



authors proposed a scaling theory suggesting that for a Coulombic {n = d — 1) two- 



dimensional diffusion-annihilation system, annihilation exponent is equal to 0.85 - close to 
but less than the mean-field exponent u = 1. Oshanin et al. suggested that the expo- 



nent should be exactly 1; Ispolatov and Krapivsky [^] came to the same conclusion using an 



" inpenetrable domain" scaling theory. In their theoretical and computational study of the 
defect annihilation in two-dimensional XY-model, Yurke et al. ||2^ argued that annihilation 



exponent is 1 , with logarithmic correction due to the logarithmic dependence of the mobility 
on the defect size. Thus, the problem of annihilation behavior for Coulombic system in two 
dimensions has not been completely resolved, although it seems rather plausible that the 
final asymptotics is governed by a classical exponent 1, with possible logarithmic correction. 
It seems to be certain, however, that annihilation in three-dimensional Coulombic system 
has a mean-field type final asymptotics with exponent 1. 

Finally, we should point out that until recently, there has been no discussion of LRI other 



than Coulombic. Recently, several scaling theories have been suggested to analyze arbitrary 
power-law interaction (more short-ranged than Coulombic) p3| , |25| , p7| .Such problems may 
arise, e.g., in describing interactions between vacancies and interstitials in a two- or three- 
dimensional crystal. 

In this paper, we propose a self-consistent theory based on coarse-grained hydrodynamic 
equations of motion for particle number density and charge density fields. This theory 
allows us to systematically calculate the dependence of the density decay law on initial 
conditions and to investigate the role of LRI. It can be shown ||3l| that the self-consistent 
approximation that we employ here is equivalent to a resummation of an infinite class of 
Feynman graphs which take into account the conserved charge density fluctuations but ignore 
the less important nonconserved number density fluctuations. This approximation can be 
further systematically improved by the use of perturbation theory and the renormalization 
group analysis that will be a subject of future publication. In the limit of weak long- 
range interactions, this approach agrees well with the known theoretical results (see Refs. 
P JT^J25| -pD|) for the two-species annihilation A + B — > 0, thereby further clarifying the 
underlying assumptions that led to these results. 

The paper is organized as follows: in section ^ we define all variables and present the 



equations of motion on which our analysis and results will be based. In section |ITl| the 
self-consistent approximation is described and the asymptotic solutions are obtained and 



analyzed. Finally, in section |I^ we analyze the resulting phase diagram and discuss it in the 
context of the previously obtained results for various diffusion-annihilation systems, which 
are selected points on our phase diagram. 



II. EQUATIONS OF MOTION 

Let us consider a system consisting of two kinds of particles, A and B, with A having a 
positive charge +q and B having a negative charge —q. We label their, time- and position- 
dependent concentrations as ni(r, t) and n2(r, t), respectively, and impose the condition that 
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< ni(r, t = 0) > = < n2(r, t = 0) > = no- The equation of motion for the densities is based 
on the generahzed law of mass conservation, violated by the annihilation process, 

dni{r,t) 



dt 

where the mass current is given by 



V- Ji = -/Cni(r,t)n2(r,t), (5) 



Ji = -DVuiir, t) - fiqMr, t)VV{r, t), (6) 

and 

is the electrostatic long-range potential at a point r at time t due to local charge fluctuations 
away from neutrality; n is the power exponent of the long-range force, /i is particle mobility 
taken to be a constant, and q is particle charge. Equations - (0) should be solved in 
conjunction with initial conditions for ni{r,t = 0) and n2{r,t = 0). Since here we are 
interested in statistical averages, rather than in a dynamic solution for a given system with 
specific initial conditions, we will focus on the density correlation functions, with averages 
over random t = initial conditions. 

Equations @ - (0) represent the coarse-grained continuum limit for the "real" equations 
of motion; all variables in these equations are averaged over "elementary volume" Uq^. Thus, 
only long-wavelength modes are actually described by these equations, and, therefore, only 
intermediate-time and large-time regimes can be analyzed. 

It is more natural and convenient to describe the system in terms of the particle number 
and charge densities. We denote the former one as p(r, t), and the latter one as f{r,t) and 
relate them to densities rii and n2 as follows: 

p(r,t) = i(ni(r,t)+n2(r,t)) , (8) 
/(r,t) = ^K(r,t)-n2(r,t)). . (9) 



Rewriting equations ® - (0) using densities / and p, we obtain: 

- DV'pir, t) = -}C{p\r, t) - f^r, t)) - gV(/(r, t) V / , (10) 

- DV'fir, t) = -gV(p(r, t) V / d'r'^I^^^ . (11) 
where Q = ^fj,q^. 

In the absence of LRI, equations ( p!0|) - ( pT]) are those analyzed in Ref. |l5|, [Tll , p!8| , |30 . 
However, the presence of the additional long-range interactions, makes the analysis of their 
asymptotic solutions rather nontrivial and, as we will show below, leads to new dynamic 
regimes. 



It is important to note that equations (|T0]) - ([TTD do not contain noise terms in their 
right-hand sides. It has been rigorously shown that such noise terms represent important 
correlations and in some cases may even become predominant in determining the asymp- 
totic decay rate. In equations of motion describing a near-equilibrium dynamics, powerful 
fluctuation-dissipation theorem determines the form of noise correlations. In contrast, in 
systems far from equilibrium, such as a system of annihilating particles, it can be shown 
|,PT|] that the effective hydrodynamic equations of motion derived from the fundametal 



master equations contain noise terms with very nontrivial correlations, of the form which 
could not be easily guessed a priori. Lee and Cardy [|l8l proved that in a two-species reac- 
tion without LRI, such noise leads only to the renormalization of the reaction rate /C, but 
not to the change of the scaling exponents, provided that space dimensionality d > 2. We 



have shown, in a similar fashion |^T|, that for systems with LRI, the noise has no effect 
on the asymptotic dynamics for c? > 2, if the renormalization of both /C and Q is implied. 
In addition, since equations ([T0|) - ( pT]) provide a coarse-grained description on the length 
scale larger than the interparticle spacing Pq^^^, the kinetic coefficients, e.g. /C, are effective 
coefficients that incorporate finite renormalization due to the correlations on short length 
scales. 

In order to simplify further analysis, we divide each of the equations (|IUp -([TT]) by po and 
transform everything to dimensionless variables as follows: 



P^p/pQ, f^f/Po, D ^ D(po 



It is important to notice that equation ([TT| ) is linear with respect to /, while equation ([lOD 
is quadratic with respect to / (this points to the system's invariance with respect to the 
simultaneous charge sign reversal for all particles). In the next section, we will describe the 
self-consistent approximation and its solutions. 

III. SELF-CONSISTENT APPROXIMATION AND SOLUTIONS OF 

EQUATIONS OF MOTION 

A. Self- Consistent Approximation 

We now make an important approximation in order to further simplify the analytical 



treatment of equations (p!0|)- (|Tl|) , namely, we choose to ignore the fluctuations of the particle 
density p and concentrate on the fluctuations of the conserved charge density /. This 
assumption is somewhat similar in spirit to the approach of Glotzer and Coniglio [^] for 
the problem of spinodal decomposition, and to the spherical approximation for the Ising 
model in the limit of — > oo. Unlike the "classical" mean- field approach, however, the 
proposed approximation does take into account exactly the charge density fluctuations, and 
is expected, therefore, to describe at least some of the features of the fluctuation-dominated 
kinetics. 

The justification of the proposed assumption lies in a simple observation that, while 
average particle density at any time is nonzero, so that < (p— <p>)^>/<p>^ 
is finite and likely to be small, average charge density is always zero, and therefore in 
comparison the charge fluctuations < (/— < / >)^ > are large. Thus, we expect the former 
fluctuations to be less important that the latter, and that we can approximate the particle 
number density by its average (time-dependent) value in the equations of motion without 
losing their important features. In a sense, our approximation is a generalization of an 
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argument used by Toussaint and Wilczek , in which they based their scahng decay law on 
a suggestion that < p >^ \/< >. It seems clear that the approximation of ignoring the 
number density fluctuations must breakdown at least below some upper-critical dimension 
duci since asymptotically p vanishes. In this case our approximation will be valid for d > d^c 
for all times and in systems below due it will be a good approximation up to a crossover 
time beyond which the asymptotics will be modified by the number density fluctuations. 
Systematically taking into account these additional fluctuations will be a subject of future 
work, ill 

Taking into account the above approximation, we rewrite the equation (0) in Fourier 
representation, taking p{t) as a spatially-independent but time-dependent function: 

+ De + gp(t)A;^-'^)/(k, t) = 0, (12) 

where a = d + 1 — n. 

The equation (|10|) in this self-consistent approximation is rewritten as: 

^ + /Cp^ = /C / ^ < /(k, t)/(-k, t) >, (13) 

where < . . . > denotes averaging over initial conditions. 

These equations of motion have to be supplemented with initial conditions. It is well 
known that the initial density distribution plays an important role in determining the scaling 
decay law. Although the self-consistent approximation employed here is well suited for a 
comprehensive study of the influence of initial conditions on the dynamics, here we limit our 
study to a single type of an initial condition. Throughout this paper we will focus on the 
dynamics initiated with a random Gaussian particle distribution, completely characterized 
by 

</(k,0)> = 0, (14) 
< /(ki, 0)/(k2, 0) > = A (27r)'^5('^)(ki + ks) , (15) 

which constraints the system to charge neutrality at all times and for simplicity we take 
the variance A = 1 (= Pq in physical units). For most charged experimental systems a 



more relevant initial conditions incorporate the suppression of long- wavelength fluctuations 
in charge density, which can be modeled (within e.g. the Debye-Hiickel approximation) by 
A(A;) = AoA;V(^^ + in equation (0). 

The diffusion-only (DO) case [Q = 0) and Coulombic case [n = d — 1) are the simplest 
systems with a relatively clear, yet interesting asymptotic behavior. All the intermediate 
interactions (arbitrary n and d) lead to a more complicated scaling behavior, with several 
regimes and crossovers. We will devote a subsection to each of these three cases. 



B. Systems without Long-Range Interactions 

There are two ways of approaching the limit of "no long-range interactions": by decreas- 
ing the force constant Q to or by increasing the power exponent n to inflnity (interaction 
with an effectively vanishing range). Obviously, these limits should give the same answer. 
For simplicity, we will set Q = and show that our self-consistent approximation yields the 
well-known results. [|]- [|ll[] 

p{t) - (ICt)-^ , for > 4 , (16) 
p{t) ~ {Dt)-'^/^ , for ci < 4 . (17) 

For Q = 0, the kinetic equation for / reduces to a simple diffusion equation, with the solution 

/(k,t) = /(k,0)e-^^'*. (18) 

Substituting this solution ([T8|) for f(k,t) into equation (|13D, and taking into account the 
initial condition ( p!4D - (1T5|), we obtain: 



An the exact solution expressible in terms of confluent hypergeometric functions is possible. 
p!8| It can also be easily shown that equation (0) describes the asymptotic solution of 
equation (|TP|). This is expected, since, as we argued above, the approximations made by 
Toussaint and Wilczek |^ are very similar to our self-consistent model. 
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The case of n — > oo will be analyzed in subsection |IIID| , where it will be shown that 
for all n > 1 + d/2, the decay law is asymptotically the same as for DO systems. 



C. Coulombic Systems 

In Coulombic systems, the long-range interaction is the strongest possible that one can 
achieve without making the system thermodynamically unstable (systems with interactions 
stronger than Coulombic have infinite pressure and chemical potential even if their total 
charge is zero). Because of this, one would expect the particle density decay for Coulombic 
systems to be very close or equal to the mean-field law p{t) ^ (/Ct)^^. ||33| As we describe 
below, the self-consistent approximation predicts the decay exponent u = 1 consistent with 
this expectation and with some simulations reported in the literature. pH] However, it can be 



shown [^] that correction to this mean-field like decay law can arise from the number density 
fluctuations and noise for d < 2, both of which have been neglected in the self-consistent 
theory presented here. 

Equation for the evolution of charge density / for the Coulombic systems in the self- 
consistent approximation can be exactly solved to yield: 

/(k, t) = /(k, 0) exp[-DkH - Q p(r)rfr] . (20) 

Using this solution (^0|) and the initial condition (14) - (^), p(t) can be easily shown to 
satisfy the following differential equation: 

^ + /Cp2 = /Cexp(-2Q /^(r)rfr) J -0-^eM-2DkH), (21) 

In order to find the asymptotic solutions, we introduce a new variable: 

0= fp{T)dT. (22) 
Jo 



Equation ( pT]) then transforms to: 



d'^Q .JdQV , IC 



+ /C — =exp(-2ge)- (23) 



df^ \dt j ' ' ^ '(l + 2Dt)^/2' 
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Let us find the " critical" dimension due-, above which the mean- field behavior is manifested. 



The mean-field solution for B is given by 

e(t) = iln(l + 0) + ..., (24) 

where . . . corresponds to subdominant constant terms and terms decreasing with time. By 
counting powers of t in the right-hand side and the left-hand side of equation (0), we obtain: 

power(LHS) = —2; 

power(RHS) = -f - 2|. 
Obviously, for the mean-field solution to be valid asymptotically, the power(LHS) should be 
larger than the power (RHS), which happens for systems with dimensionality larger than the 
critical dimension, 

d>d^, = 4(1 - . (25) 

If Q > /C, the asymptotic kinetics is determined by the slower process, which is the annihila- 
tion, with possibly interaction renormalized /C (implicitely assumed here), and the Coulomb 
interaction and diffusion are asymptotically irrelevant. For Q < }C, the Coulomb interaction 
has an interesting effect of continuously lowering the upper-critical dimension from 4 (for 
Q/IC = 0, in which case for d < A the diffusion dominates giving u = d/4) down to due given 
above. 

In order to analyze the kinetics when space dimensionality d is below due, we employ the 
"steady-state" approximation, which suggests that at long times the time derivative in the 



LHS of equation (pif ) is the smallest of the three terms. In this case, the equation of motion 
can be written as: 

^ = exp(-Qe)(l + 2Dt)-'^/^ , (26) 
6(0) = . (27) 

An exact solution of this equation is: 
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and 



p = e. (29) 

It can be easily shown that for large t, the asymptotic solution for the particle density decay 
(d < d^c < 4) is: 

P-^, (30) 

predicting the asymptotic decay exponent v — 1 ior Coulomb systems, as in the mean-field 
regime, although with Q- rather than /C-determined amplitude. This large-t limit is achieved 
when: 



(31) 



Q 

and it can be easily seen that in the limit of Q ^ (vanishing interactions) the transition 
time tL to this region becomes infinite, i.e. this time is never reached. 

If diffusion is faster than the deterministic Coulomb interaction-driven relaxation, i.e., 
D > Q, then for times less than t^, the annihilation is governed by the intermediate asymp- 
totics: 

SO the particle density is described by the Toussaint-Wilczek solution up to the crossover 
time tL- 

p = e ~ (Dt)-'^/^ (33) 

This intermediate asymptotics, which exists only for d < 4, reflects early times diffusion- 
dominated decay, with the slower deterministic Coulomb interaction driven classical 
decay appearing only at times later than t^- In contrast, for < Q or if d > 4, there is no 
extended intermediate regime and one should sec a quick transition to a classical decay law. 
Although within the self-consistent approximation the asymptotic 1/t decay is not affected 
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by the choice of reasonable initial conditions, the intermediate diffusion-dominated decay is 
certainly affected by our choice of random Gaussian uncorrelated initial conditions given in 
equation (p^). For instance if the screened Debye-Hiickel initial conditions are used with 
A(fc) = AoA;^/(A;^ + fc^), then this A(/c) will appear as a multiplicative kernel under the k 
integral in equation (^). It will then modify the intermediate decay exponent from d/4 
to u = {d + 2)/A (and, for d > 2, eliminating this intermediate region altogether), without 
modifying the asymptotic decay of equation (^). 

To sum up, we find that within the self-consistent approximation, Coulombic systems 
{n = d — 1) asymptotically exhibit the t^^ density decay, consistent with several scaling 
arguments and simulations. [p4| - |26|j33[ 



D. Intermediate systems 

Let us now consider the general case of long-range interactions with a power-law d—l< 
n < oo that is of shorter range (weaker) than the Coulomb interaction considered in the 
previous section. Equations (|I2|) - ( [T3| ) can be solved to yield: 

/(k, t) = /(k, 0) exp[-DkH - Qk'-" f p{T)dT] , (34) 
^ + = /C / ^ eM-'^DkH - 2Qe~'^ p{r)dr] , (35) 

where a = d — n + 1. 

Equation ( p5| ) is significantly more complicated than its analogs for either Coulombic or 
non-interacting cases. Nevertheless, it is possible to find its power-law asymptotic solutions. 
Using an asymptotic analysis analogous to that described in the previous subsection, we find 
several kinetic regimes depending on the values of d and n. These regimes depend crucially 
on the charge density relaxation mechanism, i.e., whether the LRI or diffusion determines 
the relaxation rate of /(k, t) at late times. In order to analyze the asymptotic behavior of 



the system, we again introduce the integrated density Q{t) as defined in equation (|22|) . We 



also assume a power-law for the density and, for d < 4, neglect the term dp/dt. In this case, 
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the equations of motion are: 



fjp, I- Ml. 

^ = ^ y ^ eM-'^DkH - 2Qk^-^m] , (36) 



Pit) = - . (37) 



Depending on a and c?, either the first or second term in the exponential in equation ( |36D 
dominates for large t, corresponding to either diffusive or super diffusive relaxation. We first 
assume that diffusive relaxation is prevalent and determine the conditions when it is true. 
In the case of diffusive relaxation mechanism and at large t, equation (^) can be simplified 
to yield: 



dB 

Itt 



where kd is a dimensionless constant absorbing integration over angular variables, and the 
integral becomes time-independent for large t if the second term in the exponential vanishes 
with time. If we assume the power-law dependence for p and B: 

p{t) oc , (39) 
Q{t) oc t^-^ , (40) 

then it follows from equation (|38| ) that for the predominantly diffusive system u = d/4, as 
expected. Thus, in order for this solution to be self-consistent, we must require that: 

1-z/<1-(t/2, (41) 
a> d/2, (42) 

and, from the definition of a, we determine the region where the relaxation and density 
decay are diffusion-limited: 

n>l + d/2, (43) 
d<A. (44) 
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This region is marked FD (fluctuation-dominated) in Figure 1. In it, the LRI are irrelevant 
at large t, although they may influence the density decay kinetics for intermediate t. The 
asymptotic decay law in the FD region is given by: 

p ~ (Dt)-''/^ (45) 

In the region referred to as the IR (intermediate region) and which lies below FD in Figure 1 
{d — 1 < n < l + d/2) the LRI are strong and dominate the diffusion at large t. To investigate 
the asymptotics of the decay in this region, we rewrite equation ( |5BD in the following way: 

(46) 

using O = At^~^. By solving equation (^) approximately we find two asymptotics in this 
region: 

p{t) ~ (Dt)-'^/^ (47) 

valid at intermediate times, and 

p{t) ^ (Qt)-'', (48) 

where 

d d 

1/ = = (49 

4: + d-2a2-d + 2n ^ ^ 

for asymptotically large times. 

The crossover time tc from the diffusion-dominated decay to the LRI-dominated decay 

is: 

2-d+2n 2 

D 2+d-n Q l + d/2-n. (50) 

Thus, in this region (marked IR - intermediate region - in Figure 1) the LRI accelerates 
the relaxation of the initial density fluctuations and thereby speed up the annihilation. If 
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Q < D, the diffusive relaxation and the ((i/4)-law may be observed for the intermediate t 
before transition to the super diffusive relaxation and faster decay takes place for t > t^. 

If d > 4 (region MF - mean-field - in Figure 1), spatial fluctuations become irrelevant 
and the classical kinetic-rate equation becomes asymptotically correct, so the decay law in 
this region is given by: 

Pit) ^ (51) 

as previously discussed. 

IV. SUMMARY AND CONCLUSIONS 

In the presented work we derive approximate kinetic equations for the annihilation- 
diffusion process with long-range forces. To analyze the asymptotic decay law for systems 
with (i > 2, we proposed a new self-consistent method of calculating the average particle 
density as function of time. Since the total particle density is a nonconserved order parameter 
with positive average at all times, we argued that its fluctuations are less important in 
determining dynamics of annihilation than that of a conserved order parameter - the charge 
density. This approximation self-consistently decouples two kinetic equations and make it 
possible to find the asymptotic solutions. 

In the limit of weak long-range interaction (either via taking n — > oo or Q — > 0), self- 
consistent equations of motion are reduced to those of Toussaint and Wilczek ||^. 

For Coulombic systems in more than two dimensions, our model yields the mean-field 
exponent z/ = 1, yet the role of segregation (i.e., charge density fluctuations) is important 
and cannot be simply left out. Ispolatov and Krapivsky proposed the "unpenetrable 
domain" scaling concept, which in Coulombic case results in decay exponent 1 independently 
of space dimensionality. Both their model and our self-consistent approximation neglect 
possible fluctuation modes due to the spatial variation of particle density, as well as noise, 
which can lead to some slowing down of the reaction kinetics, as indicated by simulations. 
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Elucidation of such modes and their role should require detailed account of noise and possibly 
use of renormalization group analysis when d = 2, since it is the critical dimension for 
the annihilation-diffusion problem. Since neglecting of noise in this problem appears to 



be justified for long-time asymptotics for systems with d > 2, |31[] it is possible that the 
exponent u = 1 for d > 2 Coulombic systems is exact, even though there is no experimental 
evidence to support this conclusion. 

The analysis of the self-consistent equations of motion ([T^) - ([T3[) suggests that there is 
a new region in the (n,d) phase diagram (the region labeled IR in Figure 1), in which the 
large-t asymptotics of the density decay is determined by long-range forces. The annihilation 
initally depletes positively charged region of negative particles and vice versa, and then the 
decay rate is determined by the speed of particle drift from such regions. Again, our large-t 
asymptotics here agrees with the "unpenetrable domain" theory of Ispolatov and Krapivsky, 
although their model predicts different boundaries of the IR region in the (n,d) phase space 
{n = 1 + d/2 is the upper boundary of the IR region in the self-consistent model, and n = d 
is the upper boundary of the IR region in the "unpenetrable domain model"; the lower 
boundary in both theories is the Coulombic line n = d — 1). The self-consistent theory also 
predicts a crossover from diffusion-dominated decay to the LRI-dominated decay at large 
times for the systems in this region. 

Because the self-consistent model is a semi-mean-field approximation (it completely ne- 
glects particle density fluctuations and only takes into account the concerved charge density 
fluctuations), it should be considered only as a first step. A systematic perturbative analysis 
of equations ([T0|) - (|TT]) around our self-consistent solution is needed to assess the role of left 



out number density fluctuations and noise. 



The proposed self-consistent model, its somewhat uncontrolled approximations notwith- 
standing, represents an important tool in qualitative analysis of dynamic processes in two- 
component systems with one conserved and one non-conserved variable. It predicts new 
annihilation behavior (IR regime) and new crossovers between diffusion-driven and LRI- 
driven decay regions, reproduces all known results for the annihilation problem in special 



limits, and can be used to systematically study the role of initial conditions in such processes. 
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FIGURES 



Figure 1. Phase diagram of the annihilation-diffusion reaction with long-range forces. FR 
- forbidden region (below Coulombic line) . IR - intermediate region - large-t asymptotics is 
determined by LRI. FD - fluctuation-dominated region - large-t asymptotics is determined by 
diffusion and initial fluctuations. MF - mean-fleld region - large-i asymptotics is determined 
by the kinetic-rate equation. 
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